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Abstract 

Fuel  cell  (FC)  power  plants  are  subject  to  a  number  of  possible  outage  and  derated  states  due  to  partial  or  full  failure  of  auxiliaries. 
Furthermore,  most  of  the  power  system  reliability  studies  reported  in  the  literature  assume  mean  values  of  the  particular  measure  of  re¬ 
liability.  However,  the  reliability  indices  for  grids  such  as  failure  rate,  outage  duration,  etc.  vary  for  different  time  period  due  to  weather 
conditions,  variety  of  power  demands  and  random  faults.  It  is  essential  to  obtain  the  estimation  of  reliability  under  all  environmental,  op¬ 
erational  and  loading  conditions.  This  paper  considers  above  mentioned  seasonal  variation  of  grid  reliability  indices  as  well  as  partial  or 
full  failure  of  fuel  cell  auxiliaries  for  grid-connected  PEM  fuel  cell  power  plants  (FCPPs).  In  the  paper,  a  detailed  state-space  model  of 
the  grid-connected  FCPP  is  presented  which  is  a  combination  of  proton  exchange  membrane  fuel  cell  (PEMFC)  power  plant  generation 
model  and  grid  outage  model.  The  state-space  generation  model  of  a  PEMFC  power  plant  is  formed  based  on  the  failure  modes  of  sys¬ 
tem  auxiliary  components.  As  for  the  grid  outage  state-space  modeling,  the  effects  of  weather  conditions  such  as  normal  and  adverse 
weather  are  taken  into  consideration  in  modeling  the  failure  and  repair  rates.  The  functional  relationship  between  weather  conditions  and 
transition  rates,  namely  failure  and  repair  rates  are  developed  based  on  the  fuzzy  set  theory  and  embedded  into  Markov  model  (MM). 
Simulation  results  are  obtained  for  a  5  kW  grid-connected  PEMFC  that  supplies  a  typical  residential  house  using  the  MATLAB  software 
package. 

©  2004  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

The  main  objective  of  any  utility  in  the  new  competitive 
environment  would  be  to  supply  customers  with  electrical 
energy  as  economically  as  possible  with  a  higher  degree 
of  reliability  and  quality.  Electric  utility  companies  have 
been  making  every  effort  to  achieve  this  objective  in 
many  ways,  one  of  which  is  to  widen  distributed  gen¬ 
eration  (DG)  usage.  Among  the  various  types  of  DGs, 
fuel  cells  (FC),  particularly  proton  exchange  fuel  cells 
(PEMFCs),  generated  tremendous  interest  for  electricity 
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and  heat  generation  due  to  its  low  operating  tempera¬ 
ture,  fast  start  up  characteristics  and  ecological  cleanness 
[1]. 

Studies  that  quantify  power  system  reliability  have  so 
far  been  limited  to  constant  transmission  rates,  covering 
two  weather  conditions,  namely  normal  and  adverse  weather 
[2,3].  In  reality,  however,  the  probability  of  a  system  or  a 
component  failure  varies  from  time  to  time  dependent  upon 
factors  such  as  a  change  of  environmental  conditions,  de¬ 
mand  variation  and  random  failures  in  the  system.  For  better 
system  operation  and  to  aid  the  formulation  of  investment 
policy,  the  variation  of  reliability  under  different  weather 
conditions  should  be  known.  Identifying  the  system  relia¬ 
bility  variations  for  different  conditions  is  crucial  for  fuel 
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cell  systems  since  it  is  an  emerging  area  of  technology.  To 
increase  competitiveness  and  market  value  of  the  FCPP,  it 
is  important  to  analyze  FC  System  Reliability  for  possi¬ 
ble  environmental,  operational  and  loading  conditions.  Al¬ 
though  numerous  studies  have  been  performed  in  the  field 
of  power  system  reliability,  to  the  best  of  the  authors  knowl¬ 
edge,  no  work  has  been  reported  in  the  area  of  FC  system 
reliability,  possibly  due  to  the  early  stage  of  FC  technol¬ 
ogy  development  and  unavailability  of  sufficient  data.  This 
paper  develops  a  grid-connected  FC  system  reliability  an¬ 
alytical  model  based  on  MM  for  multiple  weather  condi¬ 
tions  of  the  grid  and  state  of  health  of  the  FC  system.  The 
proposed  model  includes  a  detailed  state-space  generation 
model  for  FCPP  and  a  state-space  Markov  model  (MM)  for 
grid  outage.  Fuzzy  set  theory  is  used  in  the  MM  that  is  called 
fuzzy  Markov  model  (FMM),  to  incorporate  both  transition 
rates  and  temperature-based  seasons,  which  include  multi¬ 
ple  weather  conditions  such  as  normal,  less  stormy  and  very 
stormy. 

2.  Overview 


has  auxiliary  components  such  as  water  pumps,  cooling  fans, 
boiler  and  turbine.  Each  of  the  auxiliary  failure  may  produce 
many  partial  outages,  which  results  in  either  derated  or  outage 
states. 

The  number  of  outage  state  is  given  by  2n  —  1  in  Markov 
processes,  where  n  is  the  number  of  elements  that  affect 
the  system’s  state  of  health.  However,  the  probability  of 
more  than  one  auxiliary  component  failure  at  the  same 
time  is  extremely  rare  and  therefore  such  states  can  be 
neglected. 

The  general  state-space  model  of  a  generating  unit  for 
the  aforementioned  case  is  given  in  Fig.  1,  where  X  rep¬ 
resents  failure  rate,  fi  represents  repair  rate,  indices  i=  1, 
2,  . . .,  k  represent  the  number  of  components  that  cause 
unit  failure,  and  j=  1,  2,  . . .,  m  represent  the  number  of 
components  that  affect  derating  level  of  the  unit.  In  var¬ 
ious  cases,  some  other  factors  such  as  climate  conditions 
and  poor-quality  fuel  may  cause  derated  states  for  gener¬ 
ating  units.  For  instance,  fuel  quality  is  vitally  important 
for  FCPP  in  order  to  operate  them  economically  and  effi¬ 
ciently. 


2.7.  Generation  unit  reliability 

Most  generation  units  requires  a  number  of  auxiliary 
equipment  and  therefore  they  are  subject  to  different  pos¬ 
sible  derated  capacity  states  based  on  the  factors  such  as 
outages  of  auxiliaries,  fuel  quality  and  environmental  con¬ 
ditions  [4].  Hence,  the  state-space  model  of  generating  unit 
is  not  limited  to  just  possible  two  conditions  of  1  and  0, 
but  can  have  different  states  such  as  operation,  derated,  fully 
faulted  or  maintenance.  For  instance,  a  thermal  power  plant 


2.2.  Continuous  Markov  model 

In  the  field  of  reliability  assessment,  Markov  processes 
use  state  transition  rates  [5].  State  transitions  in  Markov 
processes  occur  continuously  rather  than  at  discrete  time 
intervals  and  both  system  states  and  transition  rates  together 
constitute  state-space  diagram.  Hence,  Markov  processes 
can  be  easily  applied  to  power  system  reliability  since  failure 
rates  are  equivalent  to  state  transition  rates.  The  general 
structure  of  a  continuous-time,  discrete  state  equation  for  a 
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Fig.  1.  The  state-space  model  of  the  generating  unit. 
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MM  with  n  states  can  be  described  as  [6] : 
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Eq.  (1)  can  be  rewritten  as 


d 
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where  ky  represents  the  transition  rate  between  the  ith  and 
jth  states,  pft)  represents  the  probability  of  ith  state,  and  A 
represents  the  transition  matrix.  In  the  conventional  Markov 
model  (CMM),  the  transition  rates,  ky,  are  constant  valued 
and  the  states  have  two  possible  conditions — 1  and  0.  While 
the  summation  of  the  probabilities  of  states  with  logic  1  gives 
the  availability  of  the  system,  the  summation  of  states  with 
logic  0  probabilities  gives  the  unavailability  of  the  system. 
However,  the  state-space  method  is  not  limited  to  just  two 
conditions,  1  and  0,  but  components  can  have  different  states 
such  as  operation,  derated,  fully  faulted  or  maintenance  [6,7]. 


location  of  end-users.  However,  FCPPs  are  subject  to  a  num¬ 
ber  of  possible  outage  and  derated  states  due  to  partial  or  full 
failure  of  auxiliaries  [8-10].  Some  fuel  cell  failure  modes  are 
such  as  membranes  drying  out  (insufficient  humidification), 
overheating,  passages  clogging  up  with  water,  and  freezing  of 
water  in  humidification  channels.  In  this  paper,  methodology 
in  modeling  the  reliability  of  FCPP  as  well  as  grid  reliability 
model,  which  takes  into  account  the  environmental  effects 
have  been  developed,  and  individual  reliabilities  for  FCPP 
and  grid  are  calculated.  Afterwards,  the  aforementioned  two 
considerations  are  combined  using  network  representation 
techniques  to  evaluate  the  entire  system  reliability. 


3.1.  Development  of  reliability  model  for  PEMFC 


3.  Proposed  method 

As  for  any  power  system,  the  reliability  of  a  grid- 
connected  FCPP  corresponds  to  a  measure  of  the  extent  to 
which  it  supplies  the  load  requirement  with  acceptable  con¬ 
tinuity  and  quality.  To  achieve  this  goal,  power  system  relia¬ 
bility  engineers  must  incorporate  the  reliability  of  all  system 
components  under  various  environmental  and  operational 
conditions  pertaining  to  any  specified  period.  Although  en¬ 
vironmental  effects  on  the  system  reliability  should  be  taken 
into  consideration  for  grid  supply  system,  it  is  negligible  for 
commercially  available  FCPPs  because  they  are  placed  at  the 


FCPPs  basically  convert  chemical  energy  of  hydrocarbon 
fuels,  typically  such  as  natural  gas,  into  DC  form  of  electrical 
energy.  A  FCPP  mainly  consists  of  a  fuel-processing  unit 
(reformer),  FC  stack  and  power  conditioning  unit.  The  FC 
uses  hydrogen  as  input  fuel  and  produces  DC  power  at  the 
output  of  the  stack.  A  simple  representation  of  a  FCPP  is 
shown  in  Fig.  2. 

The  performance  of  the  stack  is  expressed  by  the  polar¬ 
ization  curve,  which  gives  the  relationship  between  the  stack 
terminal  voltage  and  load  current.  In  order  to  ensure  FCPP’s 
robustness  and  adaptation  to  load  variations,  the  character¬ 
istic  of  polarization  curve  should  be  kept  at  a  constant  level 
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Fig.  2.  Basic  components  of  a  fuel  cell  power  plant. 
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by  controlling  the  parameters  such  as  reactant  flow  rate,  total 
pressure,  reactant  partial  pressure,  temperature,  and  mem¬ 
brane  humidity  [11].  To  control  these  parameters,  various 
auxiliary  components  such  as  blowers  for  cooling  and  reac¬ 
tant  air,  pump  for  hydrogen  circulation  and  humidifier  for 
reactant  air  are  used  in  the  FCPP  system.  For  power  conver¬ 
sion  and  system  control,  the  power  conditioning  unit  includes 
a  DC/DC  converter,  a  DC/ AC  inverter,  an  overall  controller, 
a  transformer  and  auxiliary  energy  storage  devices  such  as 
battery  and  super-capacitor.  However,  the  number  and  types 
of  auxiliaries  may  be  different  for  different  FCPPs  depending 
on  the  applications.  In  this  research,  a  PEM  FCPP  is  consid¬ 
ered  which  includes  the  auxiliary  components  along  with  the 
input  and  output  signals  as  shown  in  Fig.  3. 

As  shown  in  Fig.  3,  the  fuel  cell  system  includes  the  fuel 
cell  stack,  plus  all  of  the  auxiliary  equipment  such  as  air 
compressors,  pumps,  humidification  equipment,  coolers  and 
control  electronics.  During  normal  operation,  fuel  cell  sys¬ 
tem  components  such  as  compressor,  fans,  pumps,  motors, 
temperature  and  humidity  sensors,  relays,  and  other  control 
electronics  could  contribute  to  a  system  failure  or  derated 
mode  for  different  reasons  such  as  ignition  of  any  leaking 
hydrogen,  material  fatigues,  wear  outs,  break  downs,  mem¬ 
brane  drying  out,  overheating,  and  freezing  of  water  in  chan¬ 
nels  [10].  If  any  of  the  components  go  beyond  its  operational 
limit,  the  corresponding  output  will  be  degraded  by  a  fac¬ 


tor.  For  instance,  insufficient  circulating  coolant  flow  due  to 
the  failure  of  coolant  water  pump  may  cause  a  reduction  in 
the  nominal.  Other  failures  may  cause  either  a  reduction  in 
nominal  output  power  or  total  system  outage. 


3.1.1.  Development  of  state-space  model  for  PEMFC 
system 

The  PEM  fuel  cell  DC  stack  voltage  including  all  irre¬ 
versibilities  may  be  expressed  as  [12] 


P's  tack  —  Topen  Fohmic  ^activation  ^concentration 


where 


Vopen  =  N0  ■  (E°  +  EX) 


a~4  RT  (ph2-^\ 

2  E  2  F  V  m2o  ) 


is  called  the  Nernst  voltage  or  reversible  voltage  that  exists 
at  no  load  condition  for  a  given  temperature  and  pressure, 


Tohmic  —  (j  T  In)  ’  ^¥C  —  Ale  ’  RfC 


is  the  resistive  voltage  loss  due  to  the  resistance  of  non-ideal 
electrodes  and  connections  and  the  resistance  to  proton  flow 
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Fig.  3.  A  PEM  fuel  cell  based  FCPP  system. 
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in  the  PEM, 


equation  [2]. 


7dc 

is  the  voltage  loss  corresponding  to  the  activation  losses  due 
to  the  rate  of  reactions  taking  place  on  the  surface  of  the 
electrodes,  and 
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is  the  voltage  loss  corresponding  to  the  voltage  change  due 
to  mass  transport  losses,  where  c  is  (. RT/2F )  for  PEMFC.  The 
parameters  used  in  the  above  equations  are  defined  as: 

No  is  the  cell  number,  Vo  the  open  cell  voltage  (V), 
E°  the  cell  emf  at  STP  (25  °C  and  1  atm)  (V),  Ag^  the 
change  in  Gibbs  free  energy,  R  the  universal  gas  constant 
(8.1 345  J  mol  - 1  K- 1 ),  T  the  temperature  of  the  fuel  cell  stack 
(K),  F  the  Faraday’s  constant  (96485  Cmol-1),  Pu2  the  hy¬ 
drogen  partial  pressure  (atm);  Pu2o  the  water  partial  pressure 
(atm),  Po2  the  oxygen  partial  pressure  (atm),  a  the  charge 
transfer  coefficient,  I^c  the  current  of  the  FC  stack  (A),  Ium 
the  limiting  current  of  FC  stack  (A),  Iq  the  exchange  current 
density  stack  (A  cm-2)  and  c  the  empirical  coefficient. 

The  total  stack  power  is  defined  as  P  =  Vstack  •  7dc- 


Vstack  —  Tstack  NV  (9) 

In  other  words,  the  reversible  open  circuit  voltage  for  a 
hydrogen  fuel  cell,  which  is  given  by  Vopen  =  (— Agf  /IF), 
will  decrease  with  increasing  temperature.  The  situation  can 
be  summarized  for  the  other  irreversibilities  in  Eq.  (3)  as 
follows: 

•  Resistive  voltage  loss,  V0hmic  tends  to  increase  at  higher 
temperature  due  to  dry  out  mode  of  the  fuel  cell  stack. 
Because  drying  of  the  membrane  leads  to  a  poor  protonic 
conductivity  of  the  PEM  membrane,  which  means  that  an 
increase  in  resistive  losses  is  more  probable  at  higher  op¬ 
erating  temperature, 

•  Activation  voltage  loss,  Vactivation  will  increase  due  to  the 
(. RT/2aF )  factor  in  Eq.  (6).  Thus  Vactivation  is  directly  pro¬ 
portional  to  temperature  T 

•  The  only  decrease  takes  place  at  concentration  voltage  loss 
due  to  the  —(RT/2F)  factor  in  Eq.  (7).  However,  the  per¬ 
centage  of  voltage  reduction  in  Vconsentration  is  smaller  than 
the  percentage  voltage  increase  in  Vactivation  due  to  the  coef¬ 
ficient  a  that  takes  a  value  between  zero  and  one  (generally 
a  =  0.5).  Hence,  a  decrease  in  the  maximum  power  level 
will  be  observed  for  higher  temperature  as  shown  in  Eq. 
(10),  where  AP  =  AV-/dc- 


3.1.2.  Cooling  system 

To  operate  fuel  cells  effectively,  a  coolant  fluid  must  be 
circulated  through  the  fuel  cell.  The  fundamental  components 
of  the  coolant  flow  loop  system  consist  of  a  coolant  pump,  a 
radiator  and  a  fan  as  shown  in  Fig.  3.  The  cooling  system  ad¬ 
justs  the  FC  temperature  at  an  efficient  operating  point.  The 
cooling  system  includes  mechanical  and  electrical  parts  such 
as  radiator,  fan,  motor  and  motor  drives.  After  a  prescribed 
period  of  operation,  the  mechanical  parts  could  experience 
deformations  such  as  cracks,  breakdowns  and  holes  in  radi¬ 
ator,  pump  and  fan  blades  due  to  fatigue  as  well  as  wear-out 
and  stresses,  which  may  result  in  insufficient  cooling.  Fur¬ 
thermore,  electrical  motor  and/or  electrical  drives  may  mal¬ 
function  due  to  failures  in  control  circuit  and  other  related 
reasons.  The  same  types  of  failures  may  also  occur  for  other 
FC  sub-systems  that  cause  a  reduction  in  the  FCPP  output. 
The  impacts  of  the  aforementioned  auxiliary  failures  are  con¬ 
sidered  herein  rather  than  the  failure  modes  to  calculate  the 
reduction  in  overall  FCPP  system  performance.  Assume  that 
there  is  insufficient  cooling  due  to  a  reduction  in  the  cooling 
system  performance,  which  changes  the  operating  tempera¬ 
ture  of  the  FCPP  by  AT,  and  the  new  operating  temperature 
can  be  expressed  as: 

T  =  T  +  AT  (8) 

The  increase  in  temperature  leads  to  a  reduction  in  the 
FCPP  output  voltage  due  to  the  Gibbs  free  energy,  which 
is  a  chemical  energy  and  converted  into  electrical  energy, 
changes  inversely  with  temperature  as  shown  in  the  following 


P=P-AP  (10) 

Consequently,  if  the  temperature  increases  due  to  insuffi¬ 
cient  cooling,  the  power  supplying  capacity  of  fuel  cell  de¬ 
creases,  which  can  be  calculated  as: 

P  -  AP  P 

RPC(%)  =  - - -  x  100  =  -  x  100  (11) 

where  RPc  is  the  percent  reduction  in  power  supplying  capac¬ 
ity  of  fuel  cell  due  to  insufficient  cooling.  If  the  FCPP  coolant 
system  completely  fails  or  if  the  failure  causes  instability  in 
the  FCPP  operation  leading  to  an  emergency  shutdown  of  the 
system,  this  phenomenon  can  be  represented  by  a  state-space 
model  as  shown  in  Fig.  4. 

3.1.3.  Humidification  system 

Appropriate  humidification  of  the  reactant  gasses  is  vi¬ 
tal  for  proper  functioning  of  a  PEM  fuel  cell  stack.  PEM 
fuel  cells  are  quite  sensitive  to  the  humidity  of  the  reactant 
gases,  and  the  excess  liquid  water  from  the  humidification 
process  that  flows  into  the  fuel  cell  can  reduce  the  output 
power  capacity  and  may  cause  stack  damage.  The  humidi¬ 
fication  system  is  expected  to  keep  relative  humidity  of  the 
reactants  at  around  100%  level  and  a  maximum  temperature 
of  60  °C,  which  is  achieved  by  injecting  water  vapor  into 
the  air  stream  of  the  humidifier.  Biichi  and  Srinivasan  [13] 
showed  that  operating  a  PEM  fuel  cell  without  humidification 
reduces  the  fuel  cell  power  supplying  capacity  by  40%.  This 
loss  of  performance  is  mainly  caused  by  internal  ohmic  re¬ 
sistance  increase  resulted  from  drying-out  of  the  membrane. 
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Fig.  4.  The  effect  of  cooling  system  failure  on  FCPP  generation  using  state-space  representation. 


The  membrane  dry-out  mode  increases  the  magnitude  of  the 
fuel  cell  stack  impedance,  expressed  as 

Rfc  =  Rfc  +  A  R  (12) 

The  reduction  level  in  FCPP  output  power  depends  upon 
the  severity  of  insufficient  humidification.  From  Eq.  (3),  it  is 
obvious  that  the  resistive  voltage  loss,  Vohmic  is  directly  pro¬ 
portional  to  FC  stack  resistance,  which  leads  to  output  volt¬ 
age  reduction  and  consequent  power  reduction.  Since  humid 
reactants  are  a  mixture  of  dry  gasses  and  vapor,  dry  partial 
pressure  results  due  to  the  difference  between  total  pressure 
and  vapor  pressure,  expressed  as 

Pj  =  Pj  ~  Pw  (13) 

Since  humidity  ratio,  0  is  expressed  as  the  ratio  of  partial 
pressure  of  water  vapor  (Pv)  to  saturation  pressure  (Ps),  then 
the  vapor  pressure  (Pv)  can  be  expressed  as 

Pv  =  0  •  Ps  (14) 

Considering  total  system  pressure  (Pj)  and  the  factors  that 
contributes  to  the  partial  pressure,  Eq.  (14)  can  be  rewritten 
as 

Ph2  =  U)  ■  pj,  po2=x-  pj,  pu2o  =  y  •  pr  (15) 

where  w,  v  and  y  are  constants  related  to  molar  masses  and 
concentrations  of  H2,  O2  and  H2O,  respectively.  Substituting 
Eq.  (14)  into  Nernst  voltage,  Eq.  (4)  rewritten  as 

n  RT  ( w  •  v/ic\  RT 
Vopen  =  E  +yy'ln( - - - )+Jfln(pT^  ^ 

From  Eqs.  (13)  and  (16),  it  is  evident  that  if  the  total  sys¬ 
tem  pressure,  Pj,  is  reduced  to  pj,  then  the  open  circuit  volt¬ 
age  and  the  corresponding  output  voltage  will  be  reduced  by 

ff  in(/?T  -  in)- 

Assume  that  RPh  represents  the  percent  reduction  in 
power  supplying  capacity  of  the  FCPP  due  to  insufficient 
humidification.  Then  the  state-space  model  of  this  event  can 
be  represented  as  shown  in  Fig.  5.  It  may  be  mentioned  that 
possible  stack  failure  due  to  lack  of  humidification  is  not 
considered  in  the  model. 


3.1.4.  Fuelling  system 

In  a  PEM  fuel  cell,  continuous  flow  of  hydrogen  must  be 
maintained  to  ensure  seamless  electrical  energy  output.  Vari¬ 
ous  types  of  hydrogen  supply  system  are  possible  depending 
on  factors  such  as  hydrogen  production,  fuel  cell  size,  and  hy¬ 
drogen  storage.  The  simplest  system  is  to  use  high-pressure 
tank,  where  the  hydrogen  flow  rate  is  controlled  through  a 
valve.  Another  possible  way  for  performing  the  hydrogen  cir¬ 
culation  is  to  use  hydrogen-circulating  pump.  The  objective 
of  the  hydrogen  flow  control  is  to  homogenize  the  pressure 
throughout  the  stack.  As  well  as  some  mechanical  wear  outs 
and  control  malfunction,  degeneration  of  pure  hydrogen  may 
be  the  cause  of  any  possible  reduction  in  state  of  health  of 
hydrogen  supply  that  leads  to  derated  states  of  FC  unit. 

For  the  fuel  cell  generation  model,  assume  that  there  is 
an  insufficient  hydrogen  supply  to  the  FC  input.  Since  Q  = 
IF  x  wj{2  =  7dc  x  t  is  the  charge  and  wj{2  is  the  amount  of 
hydrogen  in  mol,  then 

he  =  2F  x  5^1  =  2F  x  qU2  (17) 

t 

where  qj{2  =  (wu2/t),  I&c  is  the  total  stack  current,  and  qu2 
the  molar  hydrogen  usage  per  second.  If  stack  current  is  writ¬ 
ten  in  terms  of  total  stack  power,  then  I&  =P/VS tack-  Eq.  (17) 
can  thus  be  rewritten  as 

p  =  2F  •  qn2  •  Vstack  (18) 

If  there  is  a  reduction  in  maximum  designed  level  of  hy¬ 
drogen  usage  reduces  by  Aqj{2for  any  of  the  hydrogen  circu¬ 
lating  system,  then  the  corresponding  decrease  in  the  FCPP 
output  power  can  be  expressed  as 

AP  =  IF  •  A^h2  •  Vstack  (19) 

Thus,  the  percent  reduction  in  power  supplying  capacity 
of  FCPP  due  to  insufficient  hydrogen  usage  can  be  calculated 
as 

RP  (%)  =  gH;  ~  AffH;  x  1Q0  =  4^2  x  jQQ  (20) 

qn2  qn2 

The  state-space  model  shown  in  Fig.  6  illustrates  the  effect 
of  hydrogen  supply  system  failure  on  FCPP  power  genera¬ 
tion. 


Fig.  5.  The  effect  of  humidification  system  failure  on  FC  generation  in  terms  of  state-space  representation. 
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Fig.  6.  The  effect  of  hydrogen  supply  system  failure  on  FCPP  power  generation  in  terms  of  state-space  representation. 


3.1.5.  Air  circulating  system 

The  oxygen  derived  from  air  is  used  to  complete  the  re¬ 
action  in  a  FCPP.  The  air  flow  is  controlled  by  regulating  the 
input  voltage  of  the  compressor  motor  in  order  to  maintain 
desired  amount  of  oxygen.  The  molar  oxygen  usage  of  the 
FC  per  second,  qo2  for  a  given  output  power  can  be  computed 
as 

P 

qo7  =  -  (21) 

2  4FVstack 

Eq.  (15)  can  be  adopted  to  calculate  air  usage  by  multi¬ 
plying  go2with  0-21  since  21%  oxygen  is  available  in  the  air 
mixture.  However,  in  real  applications  the  blower  and  com¬ 
pressor  are  sized  such  that  the  excess  oxygen  ratio  7  is  equal 
to  2  [14],  given  by 

supplied  rate  of  O2 

y  ~  -  (22) 

reacted  rate  of  O2 

For  modeling  FCPP  system  generation,  assume  that  the  air 
supply  is  insufficient  due  to  factors  such  as  possible  partial 
failures  and  wear  outs.  Since  abundant  oxygen  is  supplied, 
FCPP  output  will  not  be  affected  unless  partial  pressure  of 
oxygen  drops  below  the  critical  level,  which  may  lead  to 
catastrophic  membrane  failure.  An  emergency  stop  function 
will  either  isolate  the  system  from  load  or  shut  down  the 
system.  Thus  possible  failure  in  air  flow  system  either  reduces 
or  leads  to  total  air  circulation  outage.  The  failure  modes 
in  air  flow  circulation  do  not  lead  to  any  derated  state  but 
results  in  FCPP  system  failure.  The  state-space  model  in 
Fig.  7  demonstrates  the  effect  of  oxygen  supply  system  failure 
on  FCPP  based  power  generation. 

3.1.6.  Energy  storage  system 

The  hydrogen  flow  rate  is  controlled  continuously  in  order 
to  follow  the  electrical  load  variations.  However,  depending 
upon  the  FCPP  system  type,  this  flow  rate  adjustment  can  be 
achieved  with  a  time  delay,  ranging  from  a  few  microseconds 
to  30  s.  Therefore,  some  type  of  energy  storage  is  needed  for 
fast  transient  response  to  meet  the  peak  load  requirement. 
Among  the  various  energy  storage  devices,  batteries  are  the 


common  choice,  and  they  have  relatively  shorter  lifetime  and 
need  periodic  maintenance.  Although  energy  storage  selec¬ 
tion  is  a  complex  process  and  needs  to  satisfy  a  number  of 
criteria,  our  objective  here  is  to  maintain  the  power  supply 
at  the  needed  level.  If  the  FCPP  output  power  is  P  and  sys¬ 
tem  peak  load  is  Ppeak,  then  based  on  80%  battery  usage,  the 
required  energy  storage  power  is  given  by 


storage 


—  1.25  X  (Ppeak  P) 


(23) 


The  energy  storage  capacity  is  measured  as 
W  =  ^storage  x  U  where  t  is  the  number  of  hours  used 
for  supplying  the  power,  Pstorage  •  To  evaluate  the  effects 
of  energy  storage  on  FCPP  reliability,  several  factors  such 
as  maintenance,  energy  storage  lifetime,  and  mean  time 
between  failures  (MTBF)  should  be  used.  If  the  total  power 
supply  capacity  of  the  system  falls  below  the  load  demand 
due  to  energy  storage  performance  loss  and/or  failures, 
the  excess  load  can  no  longer  be  supplied  and  must  be 
disconnected  from  the  system  by  using  a  smart  energy 
management  control  system.  The  percent  reduction  in  the 
FCPP  power  supply  capacity  due  to  energy  storage  failure 
may  be  calculated  as  below,  if(P  +  Tutorage)  <  demand 


Rp  =  \(P  +  Storage)  ^demandl  x 

(P  +  Storage) 


(24) 


where, 


^storage  —  ^storage  APst0rage  (25) 

denotes  the  new  battery  power  due  to  loss  of  energy  storage 
performance.  As  a  consequence,  any  possible  energy  storage 
failure  results  in  only  derated  state,  which  is  illustrated  in  the 
state-space  model  shown  in  Fig.  8. 

Other  FCPP  components  such  as  reformer,  stack,  power 
conditioner  and  transformer  are  considered  as  essential  com¬ 
ponents  for  the  system  power  supply  and  the  failure  of  any 
of  these  component  can  bring  the  system  down.  However, 
stack  performance  has  gradual  performance  deterioration 
rather  than  catastrophic  failures.  As  the  electrodes  and  elec¬ 
trolyte  become  older,  FCPP  output  voltage,  and  its  power 


Fig.  7.  The  effect  of  oxygen  supply  system  failure  on  FCPP  based  power  Fig.  8.  The  effect  of  energy  storage  failure  on  FC  generation  in  terms  of 

generation  in  terms  of  state-space  representation.  state-space  representation. 
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Table  1 

General  format  of  the  rules 

Then  Consequent  block:  Output  is  [transition 

rates  (failure  rates,  repair  rates)] 


Antecedent  block:  If  input  is  [(weather  conditions) 
and  (load  demand)  and  (...)  and  . . .] 


drops  steadily  with  time.  Since  this  is  more  important  for 
stand  alone  FCPPs,  we  did  not  consider  this  factor  in  this 
study. 

3.2.  Development  of  fuzzy  Markov  model  (FMM) 

In  this  paper,  we  performed  the  reliability  analysis  of  grid 
using  FMM,  which  is  connected  parallel  to  PEM  FCPPs. 
Weather  effects  on  failure  and  repair  rates  are  incorporated 
in  the  grid  outage  model.  In  CMM  based  reliability  analy¬ 
sis,  it  is  assumed  that  transition  rates,  Xij  are  constant  values. 
However,  for  power  system  reliability  analysis,  the  transition 
rates  varies  from  time  to  time  depending  on  various  factors 
such  as  weather  conditions  since  the  power  system  faults  and 
their  repair  time  are  related  to  such  conditions.  Thus  the  re¬ 
lationship  between  input  parameters  (which  can  be  weather 
temperature,  load  demand  and  power  quality)  and  output  pa¬ 
rameters  (such  as  failure  and  repair  rates)  can  be  mapped  by 
fuzzy  linguistic  rules.  The  generalized  structure  of  a  FMM  for 
n  states  is  depicted  in  Fig.  9,  and  Table  1  shows  its  universe 
of  discourse. 

In  Fig.  9,  both  failure  rates  and  repair  rates  are  represented 
by  fuzzy  membership  functions  [15],  which  are  derived  from 
the  system  outage  statistics.  The  different  system  states  are 
connected  by  fuzzy  rules  and  the  input  variables  to  these 
rules  can  be  weather  conditions,  or  load  demands,  or  both. 
For  example,  weather  temperature  is  considered  as  an  input 
variable,  which  represents  climate  conditions  in  fuzzy  mem¬ 
bership  functions  and  may  be  different  for  different  regions. 

The  number  of  antecedent  variables  in  the  universe  of 
discourse  can  be  increased  or  decreased  by  considering  the 
available  data  based  on  both  system  configuration  and  region 
conditions.  For  instance,  the  universe  of  discourse  given  in 
Table  1  can  be  ruled  as: 

(1)  “If  the  weather  is  stormy  then  failure  rate  and  repair  rate 
are  high,”  or 


Worse  Moderate  Good 


Fig.  9.  Generalized  structure  of  FMM  for  the  reliability  evaluation  of  grid 
load  point. 


(2)  “If  weather  is  stormy  and  temperature  is  very  low  then 
failure  rate  and  repair  rate  are  very  high”. 

Due  to  the  complexity  of  power  systems,  faults  may  not 
be  catastrophic,  and  the  operation  of  a  power  system  can  be 
graded  from  0%  to  100%  by  fuzzy  numbers  in  the  FMM.  The 
transition  rates  in  the  FMM  are  related  to  several  variables, 
such  as,  weather  conditions,  load  demand,  etc.  The  general 
representation  of  transition  rates  X//  and  the  transition  matrix, 
T  of  the  fuzzy  Markov  state-space  diagram  can  be  expressed 
as: 

hj  =  f(Wt,Tt,Lt)  (26) 
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Fig.  10.  The  general  fuzzy  logic  diagram  for  the  developed  FMM. 
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Fig.  11.  The  algorithm  for  the  proposed  FMM  solution. 


0  Xn(Wt,Tt,Lt)  •••  lu(WuTt,Lt)~ 

hi(Wt,Tt,Lt)  0  •••  hn(.Wt,Tt,Lt ) 

•  ••  •••  •••  ••• 

Kii(Wt,  Tt,Lt)  Xn2(Wt,Tt,Lt )  •••  0 

(27) 


where  t  represents  the  sampling  time,  W  represents  the 
weather  condition  (windy,  stormy,  normal,  etc.),  T represents 
the  weather  temperature,  L  represents  the  load  demand,  and 
represents  the  fuzzy  relation.  The  relationship  between 
Xij  and  Wt,  Tt,  Lt  are  connected  by  fuzzy  linguistic  rules 
as  shown  in  Fig.  10.  The  corresponding  FMM  algorithm  is 
shown  in  the  flow  chart  of  Fig.  11. 

In  Fig.  10,  once  the  rules  are  established,  a  fuzzy  logic  sys¬ 
tem  can  be  viewed  as  a  mapping  from  inputs,  Wt,  Ti ,  and  Lt  to 
outputs,  ky.  This  mapping  can  be  expressed  quantitatively  as 
Xij  =f(Wt,Tt,Lt).  The  max-min  composition  inference  method 


is  used  in  the  inference  engine  block  and  the  centroid  method 
is  used  for  defuzzification,  defined  as  [16]. 


=  fzHc<z)z  ■  d z 
/ZM c'(z)  •  dz 


(28) 


where  fic(z)  is  the  aggregated  output  membership  function 
and  Zc  is  defuzzified  system  output. 

The  solution  process  of  fuzzy  Markov  algorithm  is  as  fol¬ 
lows: 


(1)  define  input  variables  and  rules  table, 

(2)  obtain  the  fuzzy  transition  matrix  T  (Eq.  (27)), 

(3)  defuzzification  of  Step  2  matrix,  T,  given  as  Eq.  (29), 


(29) 


(4)  Constitute  state-space  differential  equation  and  obtain 
state  probability  by  solving  the  state-space  differential 
equation  for  inputs  //.  These  steps  are  repeated  for  each 
new  input  value. 


As  for  the  FCPPs,  assume  that  the  failure  rate  of  the  aux¬ 
iliary  components  increase  with  time  until  maintenance,  and 
the  failure  rate  drops  to  its  original  value  after  maintenance.  It 
is  also  assumed  that  maintenance  is  performed  at  regular  in¬ 
tervals,  and  the  repair  times  and  the  duration  of  maintenance 
increase  when  the  components  grow  older. 

The  relationship  described  in  Fig.  12  is  estimated  by 
Weibull  model,  which  is  a  practical  tool  for  modeling  compo¬ 
nent  aging  [17].  The  Weibull  distribution  and  its  failure  rate 
are  defined  as 

/CD  =  •  e~(T/a)t>  (30) 

a? 


(31) 


where  T  is  the  component  age  (or  time  from  last  overhaul), 
p  the  shape  factor  which  determines  how  the  failure  rate 
changes  with  equipment  age  (if  /3<  1  the  failure  rates  de¬ 
creases  with  age,  if  /3  =  1  the  failure  rate  is  independent  of 
age,  and  if  /3  >  1  the  failure  rate  increases  with  age),  and  a  is 
characteristic  time  interval.  Since  there  is  no  publicly  avail¬ 
able  data  associated  with  the  failures  of  FCPPs,  a  fuzzy  logic 
rule  base  system  is  formed  based  on  the  expert  knowledge  to 
determine  the  performance  loss  of  auxiliary  components  as 
illustrated  in  Table  2. 


Fig.  12.  Reliability  model  of  grid-connected  FC  system. 
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Table  2 

The  general  format  of  the  rule  for  component  state  of  health  determination 


273 


Antecedent  block:  If  input  is  [(component  age) 
and  (maintenance  cycle)  and  (...)] 


Then  Consequent  block:  Output  is  (state  of  health  of 

component/degree  of  failure  severity) 


Once  the  individual  FC  and  grid  load  point  reliabilities 
are  obtained,  the  results  are  combined  using  network  repre¬ 
sentation  technique  for  different  operating  time  period  with 
respect  to  weather  conditions. 

3.3.  Reliability  calculation 

Normally,  MM  based  system  reliability  is  calculated  by 
summing  up  all  the  operating  state  probabilities  as  shown 
below: 

n 

i= 1 

where  Pf  is  a  row  vector  and  shows  operating  state  probabil¬ 
ities.  However,  our  concern  is  related  to  both  up  and  derated 
states.  If  these  state  probabilities  are  grouped  in  the  vector, 
then 


correction  vector,  Q,  defined  as 

Ci  =  [  1  1  •••  1  Cl  c2  ■■■  cn  flxn  (34) 

where,  c/  =  1  —  RP?%  and  RP?%  is  the  percent  power  reduction 
corresponding  to  each  derated  state. 

Hence,  the  individual  reliability  of  the  system  elements 
can  be  calculated  as 

R  =  [Piilxn  '  [Qlnx  1  (35) 

Afterward,  individual  FCPP  and  grid  reliabilities  are  com¬ 
bined  to  obtain  the  whole  system  reliability  using  network 
representation  technique.  Since  this  paper  deals  with  grid- 
connected  PEM  FCPP  that  supplies  a  typical  residential 
house  through  a  transformer,  the  reliability  model  of  the  sys¬ 
tem  can  be  depicted  as  shown  in  Fig.  12. 

Combining  the  series  and  parallel  connections,  the  system 
reliability  can  be  calculated  as 

Rsys  =  [1  -  (1  -  RGnd)  *  (1  -  RFC)]  ■  ^Tr  (36) 


Pi  =  [P?  PI 


pu  pd 

r  m  rm+ 1 


pd 

1  m+2 


pd  i 

n  J  i 


xn 


(33) 


Consequently,  the  system  reliability  analysis  throughout 
the  operating  time  can  be  made  by  following  the  methodology 
repeatedly  for  multiple- weather  conditions. 


where  “u”  represent  the  up  states  and  “d”  represent  the  der¬ 
ated  states.  Derated  state  implies  that  the  unit  may  not  be  able 
to  develop  full  rated  capacity.  Thus,  the  derated  state  proba¬ 
bilities  must  be  reduced  by  a  reduction  factor  appropriate  to 
the  deratings.  To  take  this  effect  into  calculation,  consider  a 


4.  Example 

In  this  study,  an  approach  for  reliability  analysis  of  grid- 
connected  PEM  FCPPs  is  proposed.  For  illustration  purposes, 


BULK  POWER 
SYSTEM 


Line  1 


Line  2 


lOkVA 
120V /240V 


LOAD 

(RESIDENTIAL 

HOUSE) 


Fig.  13.  Considered  grid-connected  PEM  fuel  cell  system  for  reliability  evaluation. 
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Table  3 

A  5  kW  PEM  FCPP  parameters 


Parameter 

Value 

Parameter 

Value 

Cell  Number 

88 

Oxygen  partial  pressure  (atm) 

2.17 

Cell  emf  at  STP  (V) 

0.6 

Charge  transfer  constant 

0.5 

Universal  gas  constant  (J  mol- 1  K- 1 ) 

8.1345 

Stack  current  (A) 

94.69 

Stack  temperature  (K) 

353 

Fimiting  stack  current  (A) 

105 

Faraday’s  constant  (C  mol-1) 

96485 

Exchange  current  density  (A  cm-2) 

10  6.912 

Hydr.  partial  pressure  (atm) 

1.087 

Cone,  voltage  coefficient 

0.0147 

Water  partial  pressure  (atm) 

0.464 

Stack  internal  resistance  (Ohm) 

0.00303 

Table  4 


Rule  table  for  the  functional  relationship  between  input  and  output 
parameters 


Antecedent  block 

Consequent  block  (%  reductions 
in  auxiliary  performance) 

If  age  is  Ai  (very  young) 

Then 

Reduction  is  Bi  (very  healthy) 

If  age  is  A2  (young) 

Then 

Reduction  is  B2  (healthy) 

If  age  is  A3  (medium) 

Then 

Reduction  is  B3  (medium) 

If  age  is  A4  (old) 

Then 

Reduction  is  B4  (worn) 

If  age  is  A5  (very  old) 

Then 

Reduction  is  B5  (very  worn) 

the  proposed  approach  is  applied  to  the  system  shown  in 
Fig.  13. 

A  5  kW  Plug  Power  PEM  FCPP  installed  at  the  Univer¬ 
sity  of  South  Alabama  is  considered  in  the  sample  system 
and  the  PEMFC  parameters  are  given  in  Table  3.  The  pro¬ 
posed  FCPP  generation  system  model  incorporates  the  input 
values  (component  age  and  maintenance  cycle)  and  deter¬ 
mines  the  degree  of  failure  and  failure  rates  via  membership 
functions  and  Weibull  distribution  function.  It  is  assumed  that 
the  possible  wear  out  period  starts  after  5000  h  of  operation. 
The  maintenance  process  is  repeated  at  regular  maintenance 
intervals.  However,  the  repair  time  and  the  duration  of  main¬ 
tenance  increase  at  each  cycle  as  the  auxiliary  components 
become  older. 

The  input  and  output  membership  functions  for  an  FCPP 
generation  system  modeling  are  depicted  by  the  triangular 
membership  function  shown  in  Fig.  14.  Since  maintenance  is 
done  at  regular  intervals,  component  age  gives  direct  informa¬ 
tion  about  the  number  of  maintenance  performed.  Therefore, 
the  input  membership  function  is  chosen  only  for  component 
aging. 

The  rules  used  in  the  fuzzy  inference  system  are  listed 
in  Table  4,  which  shows  the  functional  relationship  between 


Fig.  15.  The  state-space  model  of  a  5kW  PEMFC  generating  unit 
for  5th  year  operation,  where  =  1.1423  x  10-4,  A.2  =  1.1422  x  10-4, 
Z3  =  1.1421  x  1(T4,  Z4  =  1.1419  x  10-4  (failureh-1)  and  /xi=0.3088, 
fi 2  =0.3388,  /x3  =  0.3758,  /x4  =  0.4  (repair  h-1). 

component  age  and  the  reductions  in  auxiliary  performance 
level. 

Fig.  15  shows  the  state-space  model  of  a  5  kW  PEM  FCPP 
generation  unit  assuming  5  years  of  operation  is  established 
such  that  non-healthy  states  take  place  after  a  typical  5000- 
h  operation  time  and  maintenance  is  performed  at  1-year 
regular  intervals.  After  maintenance,  it  is  assumed  that  the 
components  become  similar  to  a  new  one.  In  this  model,  bat¬ 
tery  backup  is  not  taken  into  account  since  most  grid-parallel 
PEM-FCPPs  do  not  include  any  backup  energy  storage. 

In  Fig.  15,  the  performance  loss  for  oxygen  supply  system 
is  merged  with  full  capacity  available  state  since  any  reduc- 


(b) 


5  10  15  20  25 

Auxiliary  State  of  Health  Reduction  [%] 


Fig.  14.  Input  (a)  and  output  membership  (b)  functions  for  fuel  cell  generation  system  modeling. 
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Fig.  16.  Fuzzy  Markov  state-space  diagram  for  the  grid  given  in  Fig.  13. 


tion  in  oxygen  supply  does  not  lead  to  any  derated/down  state 
in  a  5 -year  operation  cycle.  The  transition  probabilities  in  the 
model  are  estimated  for  each  sub-system  based  on  the  Weibull 
distribution. 

Fig.  16  shows  eleven  operating  modes  related  to  the  state 
of  load  bus  service,  and  outage  data  for  a  1-year  basis  are 
listed  in  Table  5. 

The  input  and  output  fuzzy  variables  of  the  FMM  in  the 
grid  reliability  study,  temperature  based  seasonal  member¬ 
ship  functions,  failure  and  repair  rates,  are  derived  by  con¬ 
sidering  the  past  data  and  expert  knowledge. 

Here,  the  input  is  weather  condition  represented  by  tem¬ 
peratures  based  on  four  seasons.  Since  spring  and  autumn  are 
very  similar,  judged  by  temperatures  and  other  environmental 


conditions,  they  are  combined.  The  outputs  are  the  transition 
rates,  namely  failure  and  repair  rates.  The  triangular  mem¬ 
bership  functions  are  selected  to  model  the  transition  rates 
for  seasonal  temperature,  failure  and  repair  rates  as  shown  in 
Fig.  17. 

The  range  of  each  output  membership  function  was  deter¬ 
mined  by  the  minimum  and  maximum  transition  rates,  while 
the  maximum  membership  point  was  determined  by  the  tran¬ 
sition  rate  considering  the  maximum  number  of  occurrence 
in  the  period.  The  membership  functions  of  failure  rates  are 
labeled  as  small  (SM),  medium  (MED)  and  big  (BIG).  For  the 
repair  rates,  the  corresponding  labels  are  little  (LT),  average 
(AVR)  and  high  (HIG)  as  shown  in  Table  6.  Some  examples 
are  shown  in  Fig.  17  (b)  and  (c)  to  illustrate  the  aforemen- 


Table  5 


The  systems  operation  modes  and  its  1-year  based  outage  data 


Operating  modes 

Outage  data 

Ref. 

Failure  of 

Permanent  outage  (failure  year-1) 

Mean  time  to  repair  MTTR  (h) 

1 .  Grid  up 

No  outage 

— 

— 

2.  Grid  up 

Tri 

1.98 

30 

3.  Grid  up 

Linei 

1.86 

35 

4.  Grid  up 

Tr2 

2.04 

35 

5.  Grid  up 

Line2 

1.86 

42 

6.  Grid  up 

Tri ,  linei 

0.588 

55 

7.  Grid  down 

Linei,  Tr2 

0.318 

67 

8.  Grid  up 

Tr2,  line2 

0.594 

65 

9.  Grid  down 

Tri ,  line2 

0.312 

68 

10.  Grid  down 

Linei,  line2 

0.582 

71 

1 1 .  Grid  down 

Tri,Tr2 

0.612 

70 
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Fig.  17.  Membership  functions  of  (a)  seasons  based  on  temperature,  (b)  failure,  and  (c)  repair  rate. 


Table  6 

The  linguistic  rule  table  of  FMM 


Weather 

conditions 

Fault  and  repair  rates 

At 

Pi 

A-2 

Pi 

. . .  Xn 

Pn 

Spr./aut. 

MEDi 

AVRi 

Summer 

SMt 

HIGi 

Winter 

BIGt 

LTt 

Spr./aut. 

med2 

avr2 

Summer 

sm2 

hig2 

Winter 

big2 

lt2 

Spr./aut. 

MED„ 

AVRW 

Summer 

SM„ 

HIG„ 

Winter 

BIG* 

lt/7 

tioned  concept.  The  range  of  transition  rates  is  scaled  to  the 
interval  [0,1]  by  using  different  scaling  factors  for  each  out¬ 
put.  The  linguistic  rules  that  connects  the  input  and  output 
variables  are  obtained  for  all  possible  transitions  as  shown  in 
Table  6. 

5.  Simulation  results 

The  simulation  studies  were  conducted  to  investigate  dif¬ 
ferent  aspects  of  the  FCPP  operation.  The  simulation  of  grid 
reliability  variations  using  FMM  is  carried  out  for  multiple 
weather  conditions  based  on  temperature,  as  shown  in  Fig.  18. 
While  instantaneous  grid  reliability  for  -15,  —5,  10,  25  and 
35  °C  weather  conditions  are  shown  respectively  from  the 
bottom  to  the  top  in  Fig.  18  (a),  reliability  variation  of  the 
grid  versus  weather  temperature  is  shown  in  Fig.  18  (b)  by 
an  increased  number  of  temperature  levels. 

For  testing  the  accuracy  of  the  FMM  approach,  the  state 
probabilities  of  the  CMM  and  the  FMM  are  compared  at 
the  15  °C  weather  conditions,  which  can  be  assumed  as  a 
common  point  for  both  CMM  and  FMM.  The  results  of  this 
test,  for  each  of  the  eleven  operating  modes,  are  given  in 
Table  7.  Although  the  accuracy  of  the  results  depends  upon 
the  accuracy  of  the  models  and  system  data,  the  correlation 
between  the  results  of  the  two  models  is  encouraging. 

The  simulation  of  FCPP  reliability  variations  using  the 
proposed  approach  is  carried  out  for  different  lifetime  of 
FCPP  operation  as  shown  in  Fig.  19.  While  instantaneous 
FCPP  reliability  for  10-year  operation  period  is  given  re¬ 


spectively  from  the  bottom  to  the  top  in  Fig.  19  (a),  reliability 
variation  of  FCPP  as  a  function  of  time  (i.e.,  year)  is  shown 
in  Fig.  19  (b). 

Fig.  20  shows  the  entire  system  reliability  variation  versus 
weather  temperature  for  different  default  value  of  transformer 
reliability  and  operational  age,  which  is  obtained  by  combin¬ 
ing  the  FCPP,  grid  and  transformer  reliabilities.  It  is  evident 
from  Fig.  20  that  the  negative  impact  of  weather  temperature 


Reliability  of  Grid 


Fig.  18.  Reliability  variation  of  grid,  (a)  instantaneous  grid  reliability  for 
—  15,  —5,  10,  25  and  35  °C  weather  conditions,  and  (b)  reliability  variation 
of  the  grid  versus  weather  temperature. 
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Table  7 


The  comparison  between  the  result  of  CMM  and  FMM  for  the  15  °C  weather 
conditions  of  FMM 


Operation  states 

The  state  probabilities 

FMM  solution 

CMM  solution 

1 

0.958764 

0.959811 

2 

0.009297 

0.009061 

3 

0.010184 

0.009928 

4 

0.011168 

0.010884 

5 

0.010250 

0.009989 

6 

0.000051 

0.000049 

7 

0.000037 

0.000036 

8 

0.000067 

0.000065 

9 

0.000034 

0.000032 

10 

0.000068 

0.000065 

11 

0.000079 

0.000074 

on  grid  reliability  is  almost  vanished  due  to  parallel  connec¬ 
tion  of  FCPP,  and  the  negative  impact  of  FCPP  aging  almost 
disappeared  due  to  grid  connection.  At  first  glance,  it  seems 
that  there  is  still  considerable  weather  temperature  and  aging 
effect  on  system  reliability.  This  is  due  to  very  short  range 


Reliability  of  PEM  Fuel  Cell 


Fig.  19.  Reliability  variation  of  FC:  (a)  instantaneous  FC  reliability  for  1,  2, 
. . .,  10  years  (from  upper  to  lower  lines  respectively),  (b)  reliability  variation 
of  the  grid  versus  year. 


(b)  Weather  Temperature  [  °C  ] 


Fig.  20.  Reliability  variation  of  system  versus  weather  temperature  for:  (a) 
transformer  reliabilities  of  0.9998,  0.99985,  0.9999  (respectively  from  bot¬ 
tom  to  top),  (b)  1,5  and  10  years  of  operation  (respectively  from  top  to 
bottom)  in  case  transformer  reliability  is  0.9999. 

of  reliability  (y  axis)  shown  in  Fig.  20  (a),  and  Fig.  20  (b), 
respectively.  Indeed,  the  aforementioned  effects  on  system 
reliability  variation  are  very  near  to  zero  (~10-5). 

If  we  consider  the  load  bus  reliability  instead  of  the  entire 
system,  the  availability  of  the  bus  would  be  0.9999856  for 
the  worst  case  scenario,  and  0.999999952  for  the  best  case 
scenario.  This  implies  that  the  unavailability  of  load  bus  will 
be  within  the  annual  interruption  time  of  5  min  to  3  s,  which 
is  an  excellent  service  availability  level. 

6.  Conclusions 

In  this  paper,  a  new  approach  for  reliability  analysis  of 
grid-connected  PEM  FCPP  is  developed  based  on  the  FMM. 
While  the  transition  rates  in  FCPP  state-space  generation 
model  is  determined  using  Weibull  distribution,  fuzzy  logic 
is  used  in  the  grid  MM  to  describe  both  transition  rates  and 
temperature-based  seasons  (continuous-time,  discrete- state 
FMM).  The  FMM  is  used  to  define  the  variation  of  grid  re¬ 
liability  for  multiple-weather  conditions  for  different  oper- 
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ational  ages.  The  reliability  of  the  entire  system  consists  of 
parallel  connection  of  grid  and  FCPP,  which  supplies  load 
through  a  transformer,  is  obtained  with  respect  to  weather 
temperature  and  operational  age  of  FCPP  The  proposed 
method  is  an  effective  tool  for  reliability  modeling  of  hy¬ 
brid  systems.  Although  this  paper  deals  with  grid-connected 
PEMFC  power  supply  system,  it  gives  a  clear  methodology 
for  reliability  modeling  of  stand  alone  FCPPs  as  well. 

Although  the  transition  rates  are  constant  in  the  CMM, 
FMM  can  handle  varying  transition  rates.  If  different  weather 
conditions  are  to  be  handled  by  a  CMM,  the  number  of  states 
in  the  model  must  increase  in  proportion  to  the  number  of 
weather  conditions.  However,  the  number  of  system  states 
in  the  proposed  FMM  is  constant  for  all  weather  conditions. 
The  result  of  multiple  weather  calculation  reveals  that  the 
major  impact  takes  place  during  low  temperature  level.  As 
for  the  aging  effects  on  PEM-FCPP,  FCPP  reliability  drops 
steadily  with  time  as  the  components  get  older.  However,  in 
grid-connected  system,  the  negative  impacts  of  the  aforemen¬ 
tioned  factors  become  negligible. 
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